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We theoretically investigate the membrane fluctuations of red blood cells with focus laid on the 
role of the cytoskeleton, viewing the system as a membrane coupled to sparse spring network. 
This model is exactly solvable and enables us to examine the coupling strength dependence of the 
membrane undulation. We find that the coupling modifies the fluctuation spectrum at wavelengths 
longer than the mesh size of the network, while leaving the fluid-like behavior of the membrane intact 
at shorter wavelengths. The fluctuation spectra can be markedly different, depending on not only 
the relative amplitude of the bilayer bending energy with respect to the cytoskeleton deformation 
energy but also the bilayer-cytoskelton coupling strength. 
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I. INTRODUCTION 

A red blood cell (RBC) has been a standard test bed 
of membrane physics due to its simple structure lack¬ 
ing a nucleus and organelles. Early studies, treating a 
RBC as a lipid bilayer bag, already rendered insight¬ 
ful understanding of various properties, to name a few, 
biconcave shapesfU |2], flickering phenomenon [3], and 
tank-treading motion [3]. Still there are certain features 
signifying the role of other substructure; a spectrin net¬ 
work (cytoskeleton) coupled to the bilayer is found to 
be responsible for the shape transformation sequence [5] , 
large-scale shape changes under shear flow [B], and cre- 
nated shapes [7]. Many attempts have been made in or¬ 
der to provide simplified descriptions of the cytokeleton 
as for its effects on non-trivial shapes and fluctuation of 

RBCs [SHU]- 

Meanwhile, theoretical works nanii analyzed the fluc¬ 
tuation spectrum measured by Zilker et al. [H] and in¬ 
troduced a picture of a composite membrane, that is, a 
lipid bilayer sparsely connected to the cytoskeleton. In 
Ref. [18| . through an empirical approach, it was claimed 
that the coupling to the cytoskeleton induces surface 
tension in such a way that the effective bending mod¬ 
ulus of the lipid bilayer undergoes an abrupt jump at 
the crossover length [13]. This claim was elaborated by 
considering the cytoskeleton explicitly as spring mesh- 
work |14] . where the authors examine the elastic energy 
of the meshwork as a function of the membrane area 
coarsely grained at the mesh size ^ and thereby conclude 
that the tension contribution vanishes suddenly at length 
scale Also, a model Hamiltonian was proposed, assum¬ 
ing infinitely strong coupling of the membrane at sites 
linked to the cytoskeleton, so as to explain the induced 
tension in the long wavelength fluctuations [Kl¬ 
in this regard it seems that a RBC as a composite 
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membrane necessarily acquires surface tension at long 
length scales. However, clear signature of the surface 
tension was not always observable, and also there ex¬ 
ists subtle variance between recent observations. Yoon 
et al. reported the shape dependence of the fluctu¬ 
ation spectrum where the tension contribution is fairly 
pronounced for spherocype-shaped RBCs but almost un- 
noticeable for discocytes. In contrast, Popescu et al. ob¬ 
served clearly the induced surface tension irrespective of 
shapes |18j . Using similar interferometric technique in 
Ref. |I8j , many red blood cells showed bending dominant 
fluctuations at short wavelengths and confining dominant 
long wavelength fluctuations, but the intermediate ten¬ 
sion dominant region could be seen from few number of 
the samples [l9|. In all those experiments, the notion of 
tension jump seems to be not evident as derived from the 
theories USKIl]. The origin of such diversity is left un¬ 
known, and more experimental evidences together with 
theoretical models are anticipated in order to draw a uni¬ 
fied picture of a RBC as a composite membrane. 

In this study, we propose an exactly solvable model to 
describe elastic properties of a RBC membrane, positing 
that the cytoskeleton is spring network and the mem¬ 
brane is sparsely and regularly linked to the network via 
harmonic coupling. Unlike assumed in Ref. m, the cou¬ 
pling strength is finite and may even be weak due to 
flexibility of protein complex linking the bilayer and the 
cytoskeleton and also due to fast dissociation/association 
process of the cytoskeleton controlled by ATP concen¬ 
tration [201 Elj- Although governing factors are not 
known precisely, it is reasonable to think that the cou¬ 
pling strength can vary depending on experimental con¬ 
ditions, and observation results should intrinsically be 
diverse. In our consideration, the coupling strength is 
introduced as an essential model parameter which inter¬ 
plays with other involved energy scales to determine the 
tension emergence. 

This paper is organized as follows. In Sec. H, we 
introduce a model describing energies associated with 
fluctuations of the membrane, deformation of the spring 
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network, and the coupling between them. The effective 
Hamiltonian for the equilibrium fluctuation of the mem¬ 
brane is obtained in Sec. Ill, and analytic expression of 
the fluctuation spectrum is derived in Sec. IV. Section 
V consists of three subsections. In Sec. V. A, various 
kinds of the qualitatively distinctive fluctuation spectra 
are analyzed. We find that at length scale shorter than 
the network mesh size, the cytoskeleton elasticity is neg¬ 
ligible irrespective of the coupling strength, and the fluc¬ 
tuation spectrum of the membrane exhibits the fluid-like 
behavior. At long wavelengths, the coupling to the cy¬ 
toskeleton imparts elasticity to the membrane, and yet 
resulting features in the fluctuation spectrum thereof are 
strongly dependent on the system parameters. The effect 
of confining potential is discussed in Sec. V. B. In the 
last subsection V. C, we compare our results with previ¬ 
ous studies, in particular related to the tension “jump” 
stated in Ref. [nmi], and conclusion follows in Sec. VI. 

II. SYSTEM 

We consider a lipid bilayer of a RBC as an almost 
flat, symmetric(zero spontaneous curvature) fluid mem¬ 
brane. Assuming small deformations without overhangs, 
the membrane height can be represented by a single¬ 
valued function h{r) where r spans the flat reference 
plane at h(r) = 0. The corresponding Hamiltonian reads 
as: 

'Hm[h] — J -I-. (I) 

The second term, so called the confining potential, is in¬ 
troduced phenomenologically to describe the long wave¬ 
length fluctuations of red blood cell membranes m- 
Closed geometry [a 1 [12], periodic pinning [22] and 
excluded-volume effect HSI [23] were suggested as the 
physical origin of the confinement, yet no consensus has 
been met. Although the term has been widely used for 
many experimental studies mi ED, there are also skepti¬ 
cal views on the necessity of the confining potential mi. 
To embrace every possible situations, we do not discard 
this phenomenological term from the model. However, 
overall discussion in this paper will be made with very 
weak confinement so that the effect of the confinement is 
marginal to the fluctuation spectrum; the limit 7 —>■ O’*" 
would not hurt the key lessons of our research. More sub¬ 
tle issues on the confining potential would be discussed in 
Sec. V. B. We also mention that in order to reflect area 
conserving property of a bilayer, the bending energy usu¬ 
ally includes effective bare surface tension, which is de¬ 
pendent on the excess area of the bilayer and ambient 
temperature [211121] . However, we intentionally discard 
the bare tension in the Hamiltonian "Hm in order to ex¬ 
plicate the induced tension by the cytoskeleton. 

Another ingredient comprising the system is the cy¬ 
toskeleton which is two-dimensional regular network of 
spectrin tetramers. We consider the cytoskeleton as 


spring network |26j and model it as a sparse square ar¬ 
ray of beads connected by springs. Since we assumed the 
small deformations, the in-plane stretching and the out- 
of-plane deformation is approximately decoupled m- 
In addition, the single particle tracking experiment also 
showed that the diffusion of the band-3 proteins which 
are connected to the spectrins are significantly restricted 
Hence, we assume that the in-plane position of 
the ith bead is fixed to form a square lattice with uni¬ 
form spacing ^ and represented by the position vector 
Xi = ^(mx + ny) with integers m and n. On the other 
hand, we allow out-of-plane displacements of the beads 
and express the displacement of the Ith bead by (j){'x.i). 
The imbalance in the displacements between the near¬ 
est neighboring beads brings about either contraction or 
dilation of the connecting springs. Assuming the iden¬ 
tical spring constant k, the elastic energy cost for such 
deformation of the spring network is, 

'He[4>] = - '/'(xi + u)]*^ , (2) 

(X..u) 

where u = or u = are the primitive vectors for the 
square lattice. In fact, the cytoskeleton is closer to trian¬ 
gular (or hexagonal) network. Yet we expect that effects 
of the symmetry details must be weak on the qualitative 
properties of the membrane fluctuation and hence pursue 
simple analytic approach, remaining in the picture of the 
four-fold network. 

We finally introduce interaction between the mem¬ 
brane and the spring network. The coupling between 
the two systems is mediated by protein inclusions in the 
membrane and adaptor proteins (for example, Ankyrin) 
bridging the membrane inclusion and the spectrin net¬ 
work. Hence, the coupling must be the consequence 
of various kinds of microscopic interactions among the 
membrane, inclusions, adaptor proteins, and spectrins. 
Rather than deriving the coupling energy from all these 
complex origins, here we approximate it under the as¬ 
sumption that the separation between the membrane and 
the coupled cytoskeleton \h{xi) — (j){xi)\ is small. As the 
leading order of Taylor expansion, we can then model 
the coupling energy by harmonic potentials between the 
out-of-plane displacement of the bead and the height un¬ 
dulation of the membrane at the bead location x^: 

V[h, 4’] = Y1 > ( 3 ) 

Xi 

where v measures the coupling strength. Since recent 
numerical research based on the harmonic membrane- 
cytoskeleton interaction showed good agreement with ex¬ 
periments on red blood cell membranes [28] . we expect 
the harmonic approximation is sufficient and ignore the 
higher order terms p5| . 

In total, the system Hamiltonian of our interest is com¬ 
posed as "H = Hm +'He + V. In Ref. [H], the authors 
considered the bilayer-cytoskeleton interaction by impos¬ 
ing a hard constraint that the membrane height at the 


3 


coupling site is identical with the bead location, i.e., 
/i(xi) = (/>(xi). In our model, this corresponds to an 
extreme situation of infinite coupling (v —>■ oo), and the 
behaviors at general coupling strengths are still to be in¬ 
vestigated. Even though steric interaction can also be 
involved to correlate the membrane and the cytoskeleton 
m, this research neglects the excluded-volume effect and 
assumes the “phantom” membrane and network. 


III. EFFECTIVE HAMILTONIAN FOR 
MEMBRANE FLUCTUATION 

Equilibrium properties of the system can be extracted 
from evaluating the canonical partition function: 

Z = J V4>Vh (4) 

= J Vh , 

where the second line defines the effective Hamiltonian 
for the height deformation of the membrane, "Heft [/*■], 
which can be obtained by integrating the Boltzmann fac¬ 
tor over the fields (/'(x). 

Eor mathematical convenience, let us obtain the 
Hamiltonian in the momentum space. For the height field 
h{r) of the membrane, we employ the Fourier transform 
in two dimensions, 

= (5) 

q 

Here we discretize the space in units of the length b which 
is a microscopic length scale, e.g., molecular size of a 
lipid molecule consisting the membrane, and let the mo¬ 
mentum q take values to satisfy the periodic boundary 
conditions, = 1 for a = a;, y with Nb being the 

lateral dimension of the reference plane. We then obtain 
the bending energy of the membrane in terms of as 

Umihq] = ^E;m(q)|/iqp , (6) 

q 

where the bending energy spectrum in the discrete space 
is given by E„(q) = (4 k/& 2)[2 - cos(ga 6)]2 -f 7 . 

Since we are interested in the length scale much longer 
than the size of the monomer, we take the continuum 
limit and describe the bending energy of the bilayer as 
Emiq) = Klql"^ -by. 

The remaining components of the Hamiltonian are 
given in terms of the fields defined only at the link sites 
Xi’s, for which we introduce the Fourier transform as 

/(x,) = ^e*‘’'"‘/q/M, f = K4>. (7) 

q 

This is again the discrete Fourier transform but the 
length unit is given by ^ instead of b. In order to dis¬ 
cern the Fourier transform, Eq. ([^, we give the tilde 


symbol to the function in the momentum space. Here 
a commensurate condition ^ = pb with an integer p is 
assumed, and M = N/p is the number of link sites along 
one direction. For the periodic boundary conditions, the 
values of are given by integer multiples of 2tt jL but 
lie in the restricted Brillouin zone (RBZ), 

- <qa< tt/C , (8 ) 

where the zone boundary qc = reflects the discrete 
nature of the spring network in our theory. 

With Eq. ([^applied, the elastic energy of the spring 
network, Eq. is diagonalized into 

He[^q]= 5] i5e(q)kq|y (9) 

qGRBZ 

where X^qeRBZ denotes the summation over the range 
given in Eq. ([^, and the energy spectrum is given by 
= Ea=x,y 2 A:[l - cos(gaC]: which gives |qp be¬ 
havior only in the long-wavelength limit specified by the 
length scale qa^ 1. On the other hand, the coupling 
energy, Eq. is transformed as 

^[^q,^q]= E ^l^q-^qP- (10) 

qeRBZ 


Given the Hamiltonian components in the momentum 
space by Eqs. (§,@, and ( | 10 [ ) , we perform integration in 
Eq. Q over the field ^q and obtain the effective Hamil¬ 
tonian governing the membrane height fluctuations as 
follows: 


nes[h]=J2Em{q)\h^\^+ E 


q^RBZ 


vEejq) 

> + Ef. (q) 


\h^? • ( 11 ) 


The second term induced by the coupling to the spring 
network describes the energy for the height deformations 
occurring over length scales longer than In the inhnite 
coupling limit, one can find that the coefficient of |/iqp 
becomes i?e(q), which yields harmonic couplings in the 
real space representation, [/i(xi) — /i(xi -|- u)]^, as consid¬ 
ered in Ref. |16j . For finite coupling strength, expanding 
the coefficient as a series of 1 /u gives the first-order term, 
—and the corresponding energy is given by a 
lattice Laplacian, —[/i(xi -b u) -b li(xi — u) — 2/i(xi)]^, 
where the negative sign indicates energy gain to create 
large-scale (^) curvature. 


IV. FLUCTUATION SPECTRUM 

Now we investigate the fluctuation spectrum of the 
membrane in equilibrium, (|/iqp). The angular bracket 
denotes the equilibrium average defined as 

(|/lqP) = Z-^ j mq , 


( 12 ) 




4 


where the partition function Z and the effective Hamilto¬ 
nian Heff [/iq] are given in Eqs. Q and 0, respectively. 


In order to perform the integration, we need to express /iq 
in the effective Hamiltonian (111 in terms of /iq. Recipro¬ 
cal lattice vectors T = jQyrixX -I- Uyy) with integral 
coefficients Ux and riy are defined by the orthogonality 
relation for the discrete lattice {xi}: 


^e*'>-- = (iV/p)25q,T, (13) 


fluctuation spectra: The all-to-all coupling energy kernel 
C'(q) and the bending energy kernel for wave vectors in 
the high-order Brillouin zone i?m(q + T). Because the ef¬ 
fective cytoskeleton elastic energy CeS (q) is the harmonic 
mean of these two terms, the smaller energy contributes 
more to the fluctuation spectrum. This competition of 
two factors can result in new kinds of fluctuation spectra 
of the coupled membrane, which will be discussed in the 
next section. 


which gives the inverse transform of Eq. ([^, 


V. RESULTS AND DISCUSSION 


/iq = (M)-i^e-‘'i'"-Mx,) . (14) 

Xi 


A. Diversity of Fluctuation Spectrum 


Here we can replace h{xi) with the right hand side of 
Eq. § with r = Xi. Then, from Eq. (|T^ it follows that 


hq= P h. 


Z^'^q+T 

T 


(15) 


where q is restricted to the range in Eq. (|^. Using this 
relation, we write the effective Hamiltonian: 


E E^q+T^q(T,T')VT' , (16) 

qSRBZ T,T' 


Although the phenomenological model proposed by 
Gov et al. explains the long wavelength fluctuation data 
of Ref. [TS], the bilayer coupled to the two-dimensional 
meshwork may not follow such description for all pos¬ 
sible set of parameters k, fc, and u, and it may pos¬ 
sess potential diversity in its fluctuation spectrum. As 
given in Eq. , the characteristics of the fluctuation 
spectrum are determined by three competing energies: 
Em(q), Em(q-I-T), and C(q). Roughly speaking, for 
given q and the system parameters, one of the energy 
functions, which satisfies 


where the summation over q in the bending energy term 
is folded into the restricted Brillouin zone (RBZ), and 
the interaction kernel is given by 

Mq(T,T')=E„(q + T),5T,T'+G(q) (17) 


with 


C'(q) = 


1 vEejq) 
p^ v + Ee{q) 


(18) 


Note that C'(q) induced by the bilayer-network coupling 
describes the all-to-all coupling between any deformation 
fields with their wave vector difference given by the recip¬ 
rocal lattice vector of the spring network, reflecting the 
fact that length scales shorter than ^ are indiscernible as 
far as the network effect is concerned. 


With the effective Hamiltonian given by Eq. (161, the 
integration in Eq. (12) is performed to give 




1 - 


C(q)E-l(q) 


Am(q) ( 1 -l- C'(q) y~!T^ -|- T) 

ksT 


Am(q) + C'eff(q) 


(19) 


with 


C-ffi(q) = C-i(q) + 5] E-i(q + T) , (20) 

T^O 

which is the main results of our model for red blood cell 
membranes. Two factors enter the characteristics of the 


max[Am(q),min[C'(q),£',.,^(q-l-T)] , (21) 

determines the main properties of the fluctuation spec¬ 
trum. Here the coupling strength v regulates the form 
of C'(q). It should be also mentioned that the fluctua¬ 
tion spectrum has directional dependence which can be 
removed by angular average. We have confirmed that the 
average does not show any qualitatively different proper¬ 
ties apart from a constant factor. Every different type of 
fluctuations (e.g., non-monotonic fluctuation spectrum) 
can also be observed from the angular averaged fluctu¬ 
ation spectra. We will thus let q = {q, 0) henceforth. 

Strongly coupled membrane'. If the coupling energy v 
is sufficiently larger than the cytoskeleton elastic energy, 
i.e., u 3> 4fc = max Ee{q), the coupling energy function 
is given by 


C{q) « Ee{q)/p^ « 


( 22 ) 


where the last approximation is valid in the long wave¬ 
length limit. Figure 1 (a) shows the relevant energy func¬ 
tions for various k values. Note that in the region of in¬ 
terest, among Em{q-\-T), the term with T = —2qc makes 
the best contribution to the summation in Eq. (20). We 
only consider Exn{q—2qc) which crosses with Em(q) at the 
zone boundary q = qc = as displayed in Fig. 1(a). 
The magnitude of these bending energies in relative to 
the elastic energy depends not only on q but also on the 
spring network rigidity k, which according to Eq. (21) re¬ 


sults in the variety of fluctuation spectrum as presented 
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wave vector q 


FIG. 1: Relevant energy fnnctions and fluctuation spectra 
in the strong coupling case (w = lOOfc), where the system pa¬ 
rameters 7/(Kgc) = 10 “"^ and kl{Kq^) = 10 “^, 1 , 10 ^ are used. 
The line types for each k values are given in the panel (b). 
Competition between the energy functions presented in the 
panel (a) leads to the variety of the fluctuation spectrum as 
presented in the panel (b), depending on the network rigidity. 
The inset shows the power law exponent of the fluctuation 
spectrum, defined in Eq. (23l. The coupling to extremely 
soft network (fc = 10~^Kq^) yields insignificant modification 
of the membrane fluctuation. For k = Kq^, a smooth crossover 
from q~'^ to q~^ occurs, as the exponent in the inset displays. 
For the stiffer network {k — 10^nq^) such crossover is inter¬ 
rupted by higher-order bending energy, resulting in the non¬ 
monotonic evolution of a. 


in Fig. 1(b). We add the inset displaying the power-law 
exponent a, 


d\n{\hq\'^) 

ding 


(23) 


in order to see more clearly the wavenumber dependence 
of the fluctuation spectrum. One common feature among 
fluctuation spectra for various network rigidities is that in 
the short wavelength region (g > gc), the bending energy 
Em{q) is dominant over Em{q — 2gc) and C'(g), and the 
fluctuation spectrum shows the fluidic behavior g“^. 

On the other hand, qualitative behavior of the long 
wavelength fluctuation (g < q^) is crucially altered by the 
elasticity k. For small k/{nq^) = 10“^ (the dashed lines), 
we have that Em{q — ‘2.qc) ^ C'(g) and Em{q) > C*(g) (see 
Fig. 1(a)). Reminding of the criterion Eq. ( |2l| ), we ex¬ 
pect the fluctuation spectrum to be well described by 


the bare bending energy Em{q), which is indeed so as 
shown in Fig. 1(b). While such soft network does not 
influence the mechanical properties of the lipid bilayer, 
for large k there exists a wavenumber region satisfy¬ 
ing Em{q) < C{q) < Em{q + T) (see the curves for 
k/{nq^) = 1,10^ in Fig. 1(a)). In this case, the elastic en¬ 
ergy contribution from the cytoskeleton is observable in 
the fluctuation spectrum with q~^ behavior in Fig. 1(b). 

Note that for k/{hiq^) = 10^, near the zone boundary 
g = gc, the higher order bending energy Em{q — ‘^qc) 
is comparable or smaller than C'(g) and yields the non¬ 
monotonic fluctuation spectrum. In a recent paper [30] 
a fluid membrane discretely coupled to a very hard fluid 
membrane is found to show non-monotonic fluctuations, 
and its origin was interpreted by negative surface ten¬ 
sion. From our analysis it becomes clear that the non¬ 
monotonicity arises from the bending energy for creating 
short wavelength curvature which is less costly than the 
deformation of the extremely rigid cytoskeleton, that is, 
Em{q — ‘2qc) < C{q). The non monotonic behavior is also 
reflected on the power law exponent a which has a deep 
near q = qc, while for other k values it monotonically 
decreases from a = 4 to a = 0 with a = 2 plateau for 
k/iKql) = 1 . 

Weakly coupled membrane: Let us consider that the 
bilayer-cytoskeleton coupling is not strong enough, i.e., 
V < Ak = max i?e(q). Defining a wavenumber g„ at 
which V = Ee{qy), we obtain the approximate form of 
the coupling energy function C{q): 


C{q) 



Ee{q)/p^ , q<qv 
v/p'^ , q'^q-u 


(24) 


which well describes the behavior illustrated in Figure 
2(a). This coupling energy should be compared with 
the bending energy costs, and the fluctuation nature is 
again determined by the dominant function which follows 
Eq. (pTl). 


In the short wavelength region (g > gc), the bare bend¬ 
ing energy (g) is dominant over the other relevant en¬ 
ergy functions and results in g“^ dependent fluctuation 
spectrum, as displayed in Fig. 2(b). In the region g < gc, 
the network rigidity comes into play. For small k (see the 
dashed lines for kj{Kql) = 10“^), Ejn{q) is still a dom¬ 
inant energy factor, which explains well the fluctuation 
spectrum presented in Fig. 2(b). For rather rigid net¬ 
work, at intermediate wavenumber the coupling energy 
C'(g) can be less than Ejn{q — 2gc) and yet larger than 
Em (g). In the case, g-independent fluctuation can be ob¬ 
servable due to the coupling function behavior for g ^ g^ 
in Eq. p4| ). This coupling induced confining effect is 
more pronounced for large fc, for example, k/(nq^) = 10^. 
Such rigid network imparts the elasticity for long wave¬ 
length (g ^ qv) undulation of the membrane yielding g“^ 
behavior. This is a unique consequence of the finite cou¬ 
pling strength, which was not expected from the Gov’s 
model m and other previous studies m- 
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FIG. 2: (a) Coupling functions and (b) fluctuation spectra 
for weakly coupled membranes {v = 10~^k), where the sys¬ 
tem parameters are same as used in Fig. 1. While for soft 
networks {k/{K,q^) = 10“^, 1) the coupling does not signifi¬ 
cantly affect the fluctuation properties, extremely stiff net¬ 
work {k — lO^nq^) induces the confining effect at the inter¬ 
mediate wavelengths (q-independent spectrum) and manifest 
its elasticity in the fluctuation spectrum by q~^ dependence 
in the long wavelength region. 


B. Confining Potent iai 

Fluctuation spectrum can also be influenced by the 
confining potential in Eq. Q. Up to now, although 
we have considered only very small confining strength 7 
relative to the bending energy nq'^, the confining strength 
of red blood cell membranes is not known precisely. In 
this section, we discuss the effect of the confining poten¬ 
tial on the fluctuation spectrum. Since the confining po¬ 
tential provides a (positive) constant shift to the bending 
energy Em ( 9 ), the role of the coupling can be insignificant 
when the confinement is too strong compared to the cy- 
toskeleton elasticity. In order for the elastic contribution 
from the spring network to be visible in the spectrum, 
in particular, the q~^ dependence, the confining should 
be sufficiently weak, for which rough estimation gives a 
condition 

< VWk (for SC), 9 „ (for WC) . (25) 

Note here that \J^jk is the wavenumber at which the 
elastic energy Ef.{q) in the long wavelength limit equals 
to the confining potential, and \fkfn and q^ are roughly 


the upper boundaries that Ee(q) dominates the bending 
energies for the strong coupling (SC) and the weak cou¬ 
pling (WC) case, respectively. 

Different microscopic models have been proposed to 
understand the physical origin of the cytoskeletal con¬ 
finement. Inhomogeneous pinning due to sparse connec¬ 
tion with the cytoskeleton and the closed geometry 
of RBCs [S1IH1IT2] are two major candidates. Our results 
suggest that the two pictures entail very different physical 
situations. In our model the discrete coupling, v{h — (pY 
in Eq. ([^, indeed leads to the g-independent fluctuation 
at long wavelength. This coupling induced confinement 
is apparent only if the membrane is weakly coupled to 
very rigid cytoskeleton (see Eig. 2(b)), and in the longer 
wavelength region (9 < 9 „), one should observe the ten¬ 
sion signature as well. Eor the closed geometry of RBCs, 
the translational symmetry of the elastic network is bro¬ 
ken. This can be effectively described by adding constant 
to the elastic energy: Ee{q) —t 70 + Ee{q), which mod¬ 
ifies only the all-to-all coupling kernel C{q). Since the 
geometry induced confining 70 is now part of C{q), for 
its observability it requires the strong coupling and soft 
cytoskeleton. Such distinctive features depending on the 
origin of the confinement would help to verify which de¬ 
scription is closer to reality. 


C. Tension Emergence 


We discuss the tension emergence, comparing our re¬ 
sults with previous studies. In Refs. m 114] . the au¬ 
thors analyzed fluctuation spectra observed in Ref. m, 
and focused on the quantity K/nq where Kq, the effec¬ 
tive bending rigidity, is defined as (|/iqP) = ksT/{nqqY- 
Assuming abrupt increase of surface tension in the long 
wavelength regime, they pointed out that jump occurs 
in njKq. Our model does produce neither sudden change 
of tension nor jump in k/K q curve. Instead, the ten¬ 
sion signature only gradually appears as 9 “^-dependence 
in (|/iqp). In order to explain the tension jump in our 
model, the coupling energy function C{q) in Eq. (18) 
should vanish at 9 = 9 c. However, in our approach based 
on the microscopic model, the coupling to the spring 
network causes all-to-all correlations between membrane 


deformations, as given in Eq. (16), yielding the effective 


coupling energy C{q) which does not vanish at 9 > 90 but 
becomes less significant than the bare bending energy. 
Therefore, introducing a sudden change of tension at 9 c 
in Refs. [13 [n] cannot be justified within the physics 
explained in our model. 

Induced surface tension can also be observed in (|/iqp) 
rather than k/K q. The parameters in Refs. umi] , esti¬ 
mated from the fit to the experiment [15] . is marginal or 
discordant for observing even the gradual tension emer¬ 
gence in the fluctuation spectrum. When the bending 
modulus K = 2 X 10“^° J = 5 k^T |3 [3 [TMT5] - the 
estimated parameter values are k/{Kq^) ^ 0 ( 10 “^) and 
7 /(k 9 c) > 0(10“^), which gives > 0(1). This 
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does not meet the condition, Eq. (25), and the signal of 
induced tension in the fluctuation spectrum should be ei¬ 
ther weak or absent. When we use an order of magnitude 
larger value of the bending modulus, k = 50 which 
is also used in recent computer simulations |3T] , the non¬ 
trivial fluctuation spectrum due to the coupling with the 
cytoskeleton becomes less visible. For a given k and 7 , 
the effect of different values of k is absorbed to the dimen¬ 
sionless parameters k/nq^ and i.e., increase of k is 

equivalent to decrease of kjnq^ and ^jnq^. In Figure]^ 
the decrease of kjnq^ makes the strongly coupled mem¬ 
brane’s peculiar properties like non-monotonic fluctua¬ 
tion near q — qc and the signature of the induced tension 
be less visible. With larger bending modulus, 
gets greater so that the tension signature becomes more 
difficult to be observed. Similarly, the Figure shows 
that decrease of k/nq^ let the size of the region for wave 
vector independent fluctuation decrease. This suggests 
that, in order to observe the peculiar influences of the 
cytoskeleton on the fluctuation spectrum, the membrane 
with soft bending rigidity should be used. 

In order to observe q~‘^ behavior, k values should be 
relatively large. In fact, the spring constant can be esti¬ 
mated in various ways. From the measured shear modu¬ 
lus /r « 7 x 10“®J/m^, adopting the continuum model 
[5], we obtain the spring constant as /c = 4/i/-v/3 ~ 
10“®J/m^ [32], and correspondingly k/{Kq^) ^ 0(1). If 
treating the spectrin tetramer as an ideal entropic spring, 
the spring constant is given by fcideai = ik^T/2pLc. The 
persistence length of the spectrin p = 7.5 nm jUB inni 
and the contour length Lc = 194 ± 15 nm [34] give 
fc « 4 X 10“®J/m^ and k/{Kq^) ^ 0(10“^). These spring 
constants are orders of magnitude larger than the fitting 
values and satisfy Eq. 


brane. A membrane strongly coupled to rigid cytoskele¬ 
ton presents non-monotonic fluctuations because of the 
membrane curvature deformation at length scale smaller 
than the cytoskeleton meshsize. When soft cytoskeleton 
is coupled to the membrane, the system exhibits the well- 
studied fluctuation spectrum described by surface tension 
and the confining potential. Weakly coupled membranes 
have an extra crossover and present wave vector inde¬ 
pendent fluctuation at intermediate length scales. The 
essential source of this diversity is the competition among 
the elastic energy of the spring network, the bilayer- 
cytoskeleton coupling, and the bending energy of the 
membrane at short length scales. 

Direct measurement of such diverse spectra must be a 
difficult task, depending on instrument resolution, sam¬ 
ple availability, and control technique. Yet, as shown in 
this work, the fluctuation spectrum is dependent on the 
ratio k/nq^, and then controlling the bending modulus 
of red blood cell membranes may allow us to observe 
non-trivial spectrum like non-monotonic fluctuations for 
rigid cytoskeleton. ATP concentration can also change 
the elastic property of the cytoskeleton as well as the 
coupling strength of the anchoring proteins. More direct 
evidences should be provided by such controlled exper¬ 
iments for better understanding the physical properties 
of RBCs. Moreover, the nonequilibrium nature of RBC 
flickering in relation to active coupling between the bi¬ 
layer and the cytoskeleton is also an important issue. In 
future studies, we will pursue to analyze nonequilibrium 
spectrum in the presence of active coupling, based on our 
model with a certain type of time dependence introduced 
in the coupling strength of Eq. (|^, which, we believe, is 
a minimal approach to grasp the essential characteristic 
of the system dynamics. 


VI. CONCLUSION 

We investigate fluctuation properties of red blood cell 
membranes through the Gaussian model which takes ac¬ 
count of the finite coupling strength between the bilayer 
and the cytoskeleton. Focusing on the role of the coupling 
and discrete nature of the spectrin meshwork, we ob¬ 
tain different types of fluctuations for the coupled mem- 
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